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^ ; Abstract 

1 Bose-Einstein correlations and momentum distributions are calculated for 

r-j ', a longitudinally expanding boson source, where the expanding system has a 

finite size in space-time rapidity. Such systems are physically realized in any 

x . 

of the single jets in high energy e + e , lepton-hadron and hadron-hadron col- 
lisions as well as in high energy heavy ion collisions, where the projectile is 
not really heavy. The ID expansion generates a thermal length-scale in the 
longitudinal direction, which together with the finite length of the expanding 
tube effects both the momentum distribution and the Bose-Einstein correla- 
tion function. The Bose-Einstein correlations are shown to be more sensitive 
to the smaller, while the momentum distribution to the longer of the two 
longitudinal length-scales. 
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Introduction. There is a renewed interest in the study of relativistically expanding one 
dimensional systems partly due to the recent high energy heavy ion programme at CERN and 
AGS. The underlining space-time picture is thought to be determined by the Bose-Einstein 
correlation measurements. Both the Bose-Einstein correlations and the momentum-spectra 
were determined, see e. g. the contributions of the NA35 and NA44 collaborations in ref. [|TJ. 
In case of high energy elementary particle reactions where jets dominate the underlining 
space-time picture, usually only static parameterizations of the Bose-Einstein correlation 
functions (BECF-s) are used, with very few exceptions. For an introduction and review see 
ref. §. 

In this Letter the importance of the dynamical picture is emphasized. It seems to be 
necessary to analyze simultaneously the invariant momentum distribution (IMD) as well as 
the BECF-s in the longitudinally comoving system, LCMS, in order to measure the usually 
very long longitudinal size of the jets. 

Particle interferometry for one dimensionally expanding hydrodynamical systems were 
studied first in ref. || where Bose-Einstein interferometry was considered for an infinite, 
longitudinally expanding Bjorken-tube. Kolehmainen and Gyulassy considered the one- 
dimensionally expanding ideal inside-outside cascade ||, which was further elaborated in 
refs. f|. The correlations in between rapidity and space-time rapidity were considered 
for ID expanding strings in ref. || where also the relation of the sidewards and outward 
momentum components were clarified. However, in these works the rapidity distribution 
of the final particles remained somewhat arbitrary and there was only one characteristic 
length-scale in the longitudinal direction. This characteristic length-scale, the "length of 
homogeneity" or thermal length, appeared due to the interplay of the flow-gradient and the 
temperature in ref. || while it appeared as a characteristic rapidity - space-time rapidity 
correlation length in ref ||. In this Letter the interplay between the finite size of the boson- 
emitting source together with the finite "thermal length" shall be considered in detail. 

One-dimensionally expanding systems are thought to be realized by the hadronic strings 
in the case of 2-jet events in high energy elementary particle reactions. Any of the jets 
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of high energy pp,pp,e + e~ or other I + h collisions corresponds to an approximately one- 
dimensionally expanding system. Boson interferometry for decaying hadronic strings was 
studied first in ref. |7j and later in details in ref. ||. In these publications, however, the 
interplay between the total longitudinal length of the string and the BECF was not consid- 
ered. High energy heavy ion collisions may also create one-dimensionally expanding systems, 
especially in the case of light projectiles. Heavier projectiles may create three-dimensionally 
expanding systems which shall not be discussed here. 

Both the momentum spectra and the BECF are prescribed in the applied Wigner- 
function formalism [f|||. In this formalism the BECF is calculated from the two-body 
Wigner-function assuming chaotic particle emission. In the final expression the time- 
derivative of the Wigner function is approximated [^,[| by a classical emission function 
S(x;p), which is the probability that a boson is produced at a given x = (t, r) = (t,x,y,z) 
point in space-time with the four-momentum p = (E,p) = (E,p x ,p y ,p z ). The particle is 
on the mass shell, m? = E 2 — p 2 , and the off-shell emission functions entering into the 
BECF |J are approximated by the off-shell continuation of the on-shell emission functions. 
In this Letter the effect of final state Coulomb and Yukawa interactions shall be neglected. 

The emission function. We model the emission function in terms of the longitudinally 



boost-invariant variables: the longitudinal proper-time, r = \Jt 2 — z 2 , the space-time ra- 
pidity r] = 0.5 In 



the transverse mass m t = \ E 2 — p 2 and the momentum-space 



yt-zj 

(e + p-\ 

rapidity y = 0.5 In . We shall assume that the emission function is characterized 

by a given space-time rapidity distribution G(rj), a distribution in the transverse production 
points I{r t ) and a proper-time distribution H(t). The correlations between rapidity and 
space-time rapidity are taken into account by the conditional probability Jm t {v — y)> which 
gives the probability that a boson is emitted with rapidity y from the point-like source char- 
acterized by the space-time rapidity 77. For a thermal emission, the width of this function is 
a function of the transverse mass m t , as we shall see below. We may also introduce an arbi- 
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trary g(pt) transverse momentum distribution. Thus the emission function can be written 
in the following way: 



S(x;y,p t ) d A x = g(p t ) G{<q) H(r) I(r t ) J mt (r) - y) dr] dr d 2 r t . 



(1) 



First, we shall consider results for general, unspecified functional form of the spacetime 
distribution functions G(rj), H{t) and I{r t ) and discuss how these quantities can be deter- 
mined from the simultaneous analysis of the BECF and that of the IMD. We shall also apply 
Gaussian approximations for the distribution functions of r, 77, 77 — y and r t as follows 
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This Gaussian ansatz results in a simple analytic parameterization of the BECF as well as 
for the rapidity-dependent part of the IMD. The transverse momentum dependent factor, 
g(pt) remains unspecified and cancels from the results. 

It is to be emphasized that by introducing a finite width At] we break the boost-invariance 
of our source in the longitudinal direction too. Thus we expect a non- stationary rapidity 
distribution. The stationary rapidity distribution corresponds to the At]t — > 00 limit. 

In this model, the source at a given space-time - rapidity rj is emitting bosons with 
the rapidity y = rj, however with a non-vanishing width of the rapidity distribution of 
the emitted bosons. This width, ArjT, can be determined by Monte-Carlo simulations in 
the case of decaying strings e.g. in the Lund model At]t ~ 0.4, refs. |H||E[]. It can be 
calculated in the framework of hydrodynamics from the relativistic Boltzmann-distribution, 



p ■ u(x) 

f(x;p) = Cexp( ), where the freeze-out temperature is denoted by T and the four- 

T 

velocity of the ID expanding matter at space-time point x is given by the scaling Bjorken 
expansion [|K| as 



u(x) = (t/r, 0, 0, z/r) = (cosh(?y), 0, 0, sinh(?y)) 



(6) 



Since the four momentum p can be also written as p = (m t cosh(y),pt, m t sinh(y)), the 
relativistic Boltzmann distribution for the ID expanding systems can be rewritten as 



f{x] y, Pt) = J mt {v~y) = C exp 



m t cosh (77 — y) 



V 



C'exp 



V 



m t (y-n) 



T 



(7) 



which is the same form as in eq. [5[ Thus in the case of the thermal emission, the width of 
the rapidity distribution of the particles is a decreasing function of the transverse mass, 



Arj T (mt) 



T 



(8) 



This completes the specification of the model. 

Evaluation of the space-time integrals. We can evaluate the integrals in the LCMS, the 

longitudinally comoving system. Here we parameterize the mean momentum of the particles 

as K = (Kq, Kt,o, 0, 0) i.e. the longitudinal component of the mean momentum is vanishing. 

Since the mean momentum is purely transverse to the z axis (jet or beam axis), it is pointing 

towards a direction which is named |11,B| the out direction. This we chose to coincide with 

our x axis and index by (T,0). The y axis shall be called the side component, indexed 

with (T, S) (Transverse, Side ). In the Laboratory system the LCMS moves with the mean 
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rapidity of the pair, Y = - In 

2 



1 ' E l + E 2 +p z>1 +p z>2 \ 



V 



In the LCMS, we index the variables 



Ei + E 2 - p z ,i - p Zy2 J 

with ('). We calculate the auxiliary function S(Ak, K) in the LCMS for small values of the 
components of the relative momentum Ak = (PtQt,o,Qt,o,Qt,s,Ql)- In the integration 
we utilize the following equalities or approximations: 



rf = r] - Y and y' = y - Y (9) 
t' = t cosh(r/) « r and z' = t smh(r)') ps tt/ ~ Tory'. (10) 

The latter two approximations are valid when the width of the rapidity distribution of the 
bosons emitted at a fixed value of rj' is smaller than unity, Atjt < 1, and when the width 
of the proper-time distribution is smaller than the mean proper-time, At < r . The first 
approximation is valid for thermal emission in the m t region T < mt as well as e.g for the 
Lund strings (At/t ~ 0.4) . (Note, that in case of a decaying hadronic string, the fluctuations 
of the break-up points create the rapidity - space-time rapidity correlation length At)t-) 

The second approximation is more safe for the case of heavy ion collisions, where the 
mean freeze-out time is about r ~ 4 — 5 fm/c, and the estimated duration is very short, 
At < 2 fm/c, [|IJ. With the help of the above approximations the integrals can be evaluated. 

General results. The one-particle IMD can be calculated for arbitrary space-time dis- 
tribution functions as under the assumption that the relativistic Boltzmann distribution 
creates the correlations in between rapidity and space-time rapidity, eq. |7|. 

d n poo 

g(Pt) / drj G(r})J mt (r]-y), (11) 

J — oo 
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djdm- 

1 d 2 n 

G(y)= lim , (12) 

mt— >oo „ 

g{p t ) dydm t 

Based on the above formulas, the function G{rj) can be measured as the asymptotic (large 
m t ) limit of the rapidity dependent part of the IMD d 2 /dy/dmf. Then the measured G{rj) 
function can be convoluted with the relativistic Boltzmann distribution according to equa- 
tions 0, [TT], I2[ Thus the m t dependent width in rapidity of the IMD is prescribed from 



the model and can be contrasted to the data for a general shape of the G(rj) distribution. 
This rapidity dependent part of the IMD is obtained as a thermally broadened space-time 
rapidity distribution, and the asymptotic, large m t limit of the rapidity distribution corre- 
sponds to the space-time rapidity distribution. With decreasing m t , the thermal broadening 
becomes larger and larger except the cases when i, the rapidity plateau becomes infinitely 



long, or when ii, the decreasing part of the rapidity distribution falls outside the experimental 
acceptance. 

The general result for BECF is given in the LCMS of the particle pairs as 

C{Ak,K) = 1+ | H(PtQt,o) I 2 1 I(Qt,o,Qt,s) ? I G e Y f / Mt (r Q L ) | 2 , (13) 

Here H etc. denotes the Fourier-transformed space-time distribution and the effective space- 
time rapidity region where the boson pairs are emitted from is denoted by 

G e Y f LM = G( V )J Mt (v - Y) (14) 

dr, G( V )J Mt (v-Y) 



and the mean transverse mass is denoted by M t = (m tj i + m t)2 )/2. Thus in the one dimen- 
sionally expanding case the interpretation of the components of the correlation function is 
the following: the duration in proper-time gives contribution to the out component of the 
correlation function while the transverse distribution of the emission points, I{r t ) yields a 
symmetric contribution to both the out and the side component. The longitudinal compo- 
nent measures that region in the longitudinal direction where bosons with similar momenta 
may emerge from. 

For an arbitrary distribution of rapidities, e.g. a flat central plateau plus some wings, 
these results give an opportunity for a measurement of the underlying, non-Gaussian space- 
time distribution. 

Results for the Gaussian ansatz. The general results may be further simplified in the 
framework of the Gaussian approximations, as follows: 

d 2 n 1 y 2 

g(p t ) exp( ) (15) 



dydm 2 (2%Ay 2 (m t ) ) 1/2 2Ay 2 (m t 



and 



C(Ak, K) = l+ exp (-(R^sQ 2 T>s + R 2 t ,oQ 2 t,o + RlQl)) (16) 
where the parameters are related to those of the source by 



A V 2 + A V 2 (m t ), (17) 

R 2 g, (18) 
Rl + ^r\ (19) 

r 2 A^ 2 (m t ). (20) 

The quantity A??* corresponds to the effective correlation length between rapidity and space- 
time rapidity, as seen by the BECF. Its square is half of the harmonic mean of the squares 
of the thermal correlation length and the total size of the pseudo-rapidity distribution, 

1 1 1 

= — + , (21) 

Arjl(m t ) Ar] 2 Ar]l(m t ) 

thus it is always dominated by the shorter of the thermal and geometrical length scales 

(Ai]t and Ar], respectively). However, the square of the width of the rapidity distribution, 

Ay 2 {m t ) is the quadratic sum of the geometrical and the thermal length scales, thus it is 

dominated by the longer of the two. 

Discussion. For thermally emitting sources the Boltzmann distribution predicts a very 

specific change of the rapidity-width of the d 2 /dy/dm 2 distribution as a function of the 

transverse mass: 

T 

Ay 2 (m t ) = Arj 2 + — . (22) 
m t 

This effect, based on simple kinematics, can be experimentally checked by measuring 
the y- width of the IMD-s at different values of fixed m t . Fitting the result with the above 
equation, both the geometrical size of the emission in pseudo-rapidity, Arj and the freeze-out 
temperature T can be extracted. 

The transverse mass dependence of the width of the rapidity distribution becomes es- 
pecially well measurable in the intermediate region, where the thermal and geometrical 
length-scales are similar, see Fig. 1 and Fig. 2. as examples. 

In principle, it is possible to determine the same set of parameters from a measurement 
of the m t dependence of the Rl parameter of the BECF-s as follows: 
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This provides a possibility for the consistency check when the model is compared with 
the data. Note that this expression coincides with Sinyukov's formula for Arj — > oo, as it is 
expected since in that limit our model source corresponds to the infinite expanding Bjorken- 
tube. Thus the above expression can be considered as the generalization of Sinyukov's 
formula for finite, ID expanding systems. At present energies, we expect that the corrections 
are significant, The effects of finite geometrical size on the longitudinal radius parameter 
are visualized on Fig. 3. and Fig.4. In both cases the mean freeze-out time was chosen to be 
t = 7.0 fm/c. For Arj = 0.5 the correction at low m t yields about a factor of 0.5 decrease in 
the longitudinal radius, because the finite geometrical size is comparable with the thermal 
length-scale for moderate values of m t in the case of Fig. 3. For a larger system with 
Ar/ = 1.0 the finite size correction factor is about 0.8. The difference between Sinyukov's 
formula and its finite-size corrected version vanishes for high values of the transverse mass 
and also decreases for increasing values of the geometrical size Arj. 

The above expression is a function of three free parameters, r , T and Arj. To determine 
all the three parameters simultaneously from the m t dependence of the BECF one needs 
very dedicated set-up. It seems to be more straightforward to determine the T and At] 
parameters from the m t dependence of the rapidity - width of the d 2 n/dy/dm^ distribution, 
and use the interferometry measurement to infer the value of the mean freeze-out time To- 
Thus the importance of dedicated Bose-Einstein experiments and wide acceptance detectors 
is emphasized simultaneously: the BECF has to be studied as a function of the m t and the 
y -width of the IMD has to be determined as a function of m t , too. 

In the considered case, where the transverse position and the transverse momentum 
distribution decouple, the side component of the BECF measures the geometrical size, and 
the out component measures the quadratic sum of both the geometrical radius and the width 
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of the proper-time distribution times the mean transverse velocity of the pair. The latter 
can be considered as an increase of the effective source size in the out direction due to the 
evaporation and the mean speed of the bosons. 

Limiting Cases. Depending on the relative size of the thermal and the geometrical 
logitudinal scales, At]t and An, we have two limiting cases: Case i, If the geometrical size 
becomes much larger than the thermal length, the rapidity-width of the IMD approaches 
the geometrical size and becomes independent of the transverse mass, 

j2 2 

an y 

oc exp( ) if Ai] » An T , (24) 

dydm 2 2An 2 

T 

R\ = r 2 Ar] 2 = r 2 — if Ar] » An T . (25) 
m t 

In this limit, corresponding to long and relatively cold expanding system, the geometrical 
longitudinal size is dominating the width of the rapidity distribution, while the thermal longi- 
tudinal length-scale is dominating the longitudinal radius parameter of the BECF, measured 
in the LCMS of the boson pair. 

Case ii, If the rapidity distribution of the particles at a fixed space-time rapidity is much 
more wider than the geometrical length-scale, then interferometry measures the geometrical 
size properly, and the temperature can be inferred from the width of the rapidity distribution: 

d 2 n y 2 E 
oc exp( ) cx exp( ) if An T >> An, (26) 

dydm 2 2Ar] 2 T T 

R 2 L = r 2 An if Ar] » An T . (27) 

This limiting case corresponds to a relatively hot and small, expanding system. The 
effects of the expansion cancel from the results here. The expressions are the same as those 
for a hot source with no expansion. 



Note that this work is an application of the ideas presented in |L3[ to the case of the 
relativistic one dimensionally expanding systems. 
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In summary, we have calculated the rapidity distribution and the Bose-Einstein correla- 
tion function in the longitudinally comoving frame for one-dimensionally expanding systems. 
In the longitudinal direction, there are two length-scales present: the geometrical and the 
thermal one. The latter is either due to the combined effect of the flow-gradient and the 
temperature or e.g. to the stochastic nature of the string fragmentation. In an analytically 
tractable model for the emission function, the longitudinal component of the Bose-Einstein 
correlation function is shown to measure dominantly the shorter while the width of the 
rapidity distribution measures dominantly the longer length-scale. This confers the intu- 
itive picture that the strings are really long, but cannot be seen by BEC measurements 
due to their expansion. The geometrical size can be determined from the m t dependence 
of the rapidity-width of the one-particle invariant momentum distribution. The freeze-out 
temperature, the mean proper-time, the width of the proper-time distribution as well as 
the geometrical size in the transverse direction can be determined within this picture by 
a combined use of the radius parameters in the correlation function on one hand and by 
the analysis of the transverse mass dependent rapidity-width of the invariant momentum 
distribution on the other hand. 
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FIGURES 

Fig. 1. 




FIG. 1. The change of the rapidity-width of the cPn/dy/drrif distribution for increasing values 
of the (fixed) transverse mass, rat = 0.3,0.5,1.0 GeV/c, as shown by the solid, long-dashed and 
short-dashed lines, respectively. The dotted line stands for the asymptotic m t — > oo limit, coincid- 
ing with G(y), the geometrical shape of the source in the longitudinal direction which was chosen 
to be a Gaussian with Ai] = 0.5. The distributions were normalized to the same maximum value 
of 1. The pion mass (m = 0.14 GeV) and a freeze-out temperature of T = 0.14 GeV was utilized. 
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Fig. 2. 
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FIG. 2. Same as Fig. 1. but for At? = 1.0 . 
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Fig. 3. 
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FIG. 3. The m t dependence of the Rl radius parameter for finite ID expanding systems is 
shown by the solid curve. The mean freeze-out time is chosen to be to = 7 fm/c and At] = 0.5. 
The other parameters are the same as for Fig. 1. Dashed line shows the prediction for Rl based 
on Sinyukov's formula. 
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Fig. 4. 
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FIG. 4. Same as Fig. 3. but for Arj = 1.0 . 
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